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The pion-nucleon and pion-nucleus interactions have been much studied, both theoret- 
ically and experimentally, for many years. However, the corresponding interactions of the 
eta meson, mainly because of the lack of 77-beams, have - by comparison - been neglected. 
The main interest in 77's has been the possibility of ^-nuclear quasi-bound states. Such 
states were first predicted by Haider and Liu [jy and Li et al. ||, when it was realised 
that the ^-nucleon interaction was attractive. Calculations by Ueda indicated that this 
may happen already in the //-deuteron system. If these states exist, then one may expect 
them to be narrow in few-nucleon systems, and so be easier to detect there. The first 
verification of this hypothesis was made by Wilkin Q, who has suggested that an indirect 
effect of such a state is seen in the rapid slope of the pd — > 3 Herj amplitude detected 
just above the r\ production threshold ||. Also an indication of strong three-body rjpp 
correlations follows from the measurement of the pp — > ppr\ cross sections in the threshold 
region 0. 

It has been shown in Ref. |7[ that the strengths of rjd interactions, in particular the 
magnitude of the scattering length in this system and the position of quasi-bound or 
virtual state are very sensitive to the value of the rjN scattering length. Moreover, the 
behaviour of the 77 N scattering matrix off the energy-shell was found to be important. 
In order that these 7/-nucleus studies can be put on a firmer foundation, it is, therefore, 
necessary that a better parametrization of the basic 7/-nucleon interaction be available. 
With this in mind, a three channel analysis of the eta-nucleon(?7iV), pion-nucleon(7riV) 
and the two-pion-nucleon(7T7riV) three-body system is carried out. This is done in terms 
of a fT-matrix based on pion-nucleon amplitudes and eta-production cross sections - the 
actual data being the irN amplitudes of Arndt et al. ||, the ttN — > rjN cross sections 
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reviewed by Nefkens || and the jp — > rjp data of Krusche et al. ||10|| . In Ref. |]TTJ it is 
shown that the photoproduction cross section runs essentially parallel to the electropro- 
duction cross section in the region some 100 MeV above threshold. Therefore, including 
electroproduction data into the following analysis is not expected to lead to a different 
conclusion. 



This analysis is now carried out and, we believe, it improves analyses done directly 
in terms of resonant T matrices. The main motivation for this study is to extract the 
eta-nucleon scattering length and effective range and to determine in these quantities the 
uncertainties allowed by the existing data. The next check and refinement is expected to 
follow from the few-body rj physics. 

For s-wave scattering in a system consisting of the two channels ttN and r]N - here 
denoted simply by the indices 7r and rj - the i^-matrix, which is essentially a generalisation 
of the scattering length, and the T- matrix that follows from it, can be written as 

(1) 



K 



l K K \ 

± *- 7T7T ± *- TjTT 



and T 



1 IQn A-7T7Y 1 ifirjA-Tjrj 
Y 1 i>C[ri Arfr-j 1 iCfrjArjrj J 



where q niTj are the center-of-mass momenta of the two mesons in the two channels ir, rj. 
The channel scattering lengths Aij are expressed in terms of the .fT-matrix elements, via 
the solution of T = K + iKqT, 



Ann — K n7T + iK^q^ / [1 — iq^Kyrj), — K V7T /(1 — iq^K^ 

= K m + iK^ n q T / (1 — iq^K^). 



(2) 



These equations form a basis in which to describe two channel scattering in terms of the 
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parameters of the K-matrix. These .fT-matrices must account for several observed features 
of the experimental data - in particular: 

a) The 5-wave irN resonances 5(1535) and 5(1650). The effect of these is inserted as 
poles at E — E and E±, which are treated as free parameters. However, their values are 
expected to be near 1535MeV and 1650MeV. Differences arise, since E 0jl are renormalised 
by the presence of two background terms K nv and K m , which describe other forms of the 
interactions and channel couplings not included explicitly. 

b) Experimentally, the rj does not appear to couple to the 5(1650) resonance, and so 
this coupling is not included in the model. 

c) There is a small correction for inelasticities of the 5(1535) and 5(1650) due to cou- 
plings to the two-pion nucleon channel. This is treated in an " optical potential manner" , 
which means the introduction of the two-pion channel K matrix and its subsequent elim- 
ination. This leads to a complex correction to the two channel K matrix. For example, 
the 5(1535) has a coupling to the three-body channel described by a singular if-matrix 
if 3 3 = an d its coupling to the two-body channels by K 3<i = > where i = ir,r). 
In the three body-channel, there is a relative momentum equivalent to the above qi. This 
is a three-body phase space element q 3 . It may be included together with the coupling 
parameter 73 into a small contribution to the width of the 5(1535) width r 7r 7r /2 = 73^373. 
Only this combination of the two-pion parameters enters the correction to the K matrix 
of the two channel problem. In principle, it should be proportional to the three body 
phase space, and this energy dependence has been accounted for. Now, the correction 
to the basic two channel K matrix, which stems from the three body channel is readily 
obtained to be 
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5R o = K h3q3 K^ 
l-iq 3 K 3> 3 



A similar procedure is applied to describe the slightly higher inelasticity of the 5(1650) 
resonance. 

These features are included in the A'-matrices as follows. 



JV 7T7r Eq _ e T El _ E T ' i_ iq3K33 ) ^nr/ ^irr] ' Eq _ e T l i- iq3 K 33 ' 



wWp K - _l_ 73(1) K _ ^(0)73(0) V7. (1)73(1) 

wnere A 33 - + , - — 1 — . 



K v^ 73(0) (,) 

Kr > 3 = ~F W W 

-C/n — Hi 



In the above model, there are 10 parameters that are determined by a Minuit fit to 110 
pieces of data - 23 are ttN amplitudes (real and imaginary) ||, 11 are ttN — > r]N cross- 
sections[cr(7rr7)] || and 53 are 7iV — > r/N crossections [(7(777)] fllOl . In practice, the actual 
cross section data was used in a reduced form, from which threshold factors have been 
removed - namely: 



a(iTT]) r = a(irr))— and r{jr]) r = J a{jr])—^. (6) 
% V ^% 



The values of T{^r\) r given in Ref. [|10j are used here directly, even though the mass of 



the 77 there is 547.12MeV, compared with the present value of 547.45MeV in Ref. [[12 



Such small differences are unimportant here, since the main threshold effect is removed 



by considering the combination cr/q v . In terms of the scattering amplitudes (T) of eq. ([[]), 
the corresponding model expressions are: 

2q n 



a{jirj) r = An (Re T^ v ) + (Im T nv ) - — and r(^7]) r = A(Phot)^J (Re T m ) 2 + (Im T m f 

oq v 

where A(Phot) is a normalisation parameter that simulates the actual production am- 
plitude. This parameter is assumed to be energy independent and is treated as a free 
parameter in the Minuit minimization. The resulting fit had a \ 2 of 0.83/dof and the 
outcome is seen in Fig. 1. Since it is not clear that the four sets of data in Fig. 1 have 
equal weight, it is of interest to also look at the separate x 2 /dpt - A) 0.73 B) 0.75 C) 
0.94 and D) 0.60. This shows that, indeed, good fits are achieved in all four sets of data 
and that the overall x 2 /dof is not dominated by any particular set. 

From Table 1 it is seen that those parameters that can be compared with numbers in 



the Particle Data Tables 12 fall into three classes: 



a) T (Total), i](br), ir(br), T (Total, 1) and ir(br, 1) can be compared directly and are 
seen to be consistent with the experimental uncertainties. The relationship between the 
above T's and the 7's is determined by the T matrix, which - close to the resonance 
- should be of a Breit-Wigner form with an energy dependent width. This relates the 
channel parameters 7 to the total width T , with elasticities and the channel momenta 
calculated at the resonance energy q(PDT). Thus, for example 7^ = 0.5ir(br)T/q(PDT). 

b) Eq and E\ are the positions of the bare poles in the ii'-matrices. As mentioned 
earlier, these get slightly renormalised in going from i^-matrices to T-matrices to give the 



numbers in the Particle Data Tables 12 
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c) The seven parameters in a,b) are essentially obtained by fine-tuning the correspond- 
ing experimental numbers - as is seen by the close agreement between the two. However, 
the remaining three parameters K m , K nv and A(Phot) are completely free. In principle, 
the first two could be related to some more fundamental model based on some underlying 
lagrangian as in Refs. [fl3 1 . The third parameter could also be calculated, if a mechanism 



for r] photoproduction were used. 

The values of the branching ratios t](br), n(br) for the S"(1535) resonance also give a 
prediction for the two-pion ratio to be 1 — t](br) — 7i(br) = 0.038 - a number in line with 
experimental estimates of 0.05-0.20. 

The errors on a, r and s were obtained by repeating the calculation for a random 
selection of the nine parameters defining the .^-matrices of eq. ([!]). This selection was 
chosen to ensure the distribution of each parameter was a gaussian centered on the values 
in Table | and with the same standard deviation. Several tests were made to determine 
the dependence of these errors on the number of runs and on the size of the region each 
side of the gaussian maximum over which the random points were chosen. The errors 
shown are for 1000 runs using regions that were 3 standard deviations. 

The negative sign of the effective range is expected, since it arises quite naturally due 
to the proximity of the S^ISSS). For the single r] channel case dominated by the resonance, 
one would have r q^/2 = {EthreshoU — E)/j v . This is a fairly large negative effective range 
of about -3 fm. The presence of other channels and background terms reduce it to about 
half of this value. The shape parameter appears to be small. In fact the imaginary part 
is consistent with zero. 
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In Table 2 a comparison is made with earlier determinations of the scattering length. 
There it is seen that the present result supports, in particular, the estimates of Refs. M 



and [113]. It is difficult to compare with the other references, since they do not give any 



error estimates. 

Fig. 2 shows that, within 30MeV of the rj threshold, the effective range expansion is very 
good. For a parametrization up to lOOMeV from the threshold, the effect of the shape 
parameter(s) plays an increasingly important role. Also it is seen that the effective range 
must be included, if the rjN scattering is needed 10-20MeV away from the r\ thresh- 
old at 1485. 7MeV. Such excursions from the threshold are needed, for example, when 
extrapolating below the threshold in ^-bound state situations. In the present case, the 
threshold value of the r]N amplitude (0.75+i0.27) becomes 0.49+i0.10 fm at 1468.4MeV 
and 0.51+i0.51 fm at 1500. OMeV. Such differences could be crucial in discussions con- 
cerning the existence, or not, of r?-bound states in few nucleon systems. 

One of the authors (S.W.) wishes to acknowledge the hospitality of the Research 
Institute for Theoretical Physics, Helsinki, where part of this work was carried out. The 
authors also thank Drs. R. Arndt and B. Krusche for useful correspondence and Drs. 
J.Niskanen and M.Sainio for several discussions. 
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TABLES 

TABLE I. The optimised parameters from Minuit defining the i^-matrices and the corre- 
sponding values from the Particle Data Tables (PDT) |i~2]| . 





K 


K 


£ (MeV) 


£i(MeV) 


F(Total)(MeV) 


Minuit 


0.177(33) 


0.022(13) 


1541.0(1.6) 


1681.6(1.6) 


148.2(8.1) 


PDT 






1535(20) 


1650(30) 


150(50) 




r/(br) 


n{br) 


r {Total, l)(MeV) 


7r(6r, 1) 


A(Phot) 


Minuit 


0.568(11) 


0.394(9) 


167.9(9.4) 


0.735(11) 


19.74(36) 


PDT 


0.30-0.55 


0.35-0.50 


145-190 


0.55-0.90 
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TABLE II. Results compared with earlier works. The numbers in [...] are the values of a 
and ro, when the exact scattering amplitudes are fitted with s = 0. 



Reference 


Scattering Length(fm) 


Bhalerao and Liu [14 

Benhold and Tanabe 15] 
Arima, Shimizu and Yazaki [O] 
Svarc, Batinic and Slaus [O 
Wilkin |] 

Sauermann et al. 13 1 

.ri Uditj V a>Li\JL 1 N t-livdlo lO 


0.27+i0.22 
0.28+i0.19 
0.25+i0.16 
0.980+i0.37 
0.886+i0.274 
0.55(20)+i0.30 
0.51+i0.21 


This paper 




Scattering length(a) 
Effective range(ro) 
Shape parameter(s) 


0.751(43)+i0.274(28) 
[0.751(43)+i0.274(28)] 
-1.496(134)-i0.237(37) 
[-1.497(134)-i0.237(38)] 
-0.102(15)-i0.008(10) 
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FIGURES 



FIG. 1. The if-matrix fit to experimental data as a function of the center-of-mass energy 
Ecm: A) The ttN — ► rjN data of Ref. H - the reduced cross-section in mb containing the factor 
q-rc/qr], B) r(7^) r the reduced cross-section of Ref. [pj in units of 10~ 3 /m 7r + , C) The real part 
of the ttN amplitudes (q n Re T) H, D) The imaginary part of the irN amplitudes (q n Im T) 
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FIG. 2. The quality of the effective range expression versus the exact values. The solid line 
shows the exact results, the dashed line the effective range expansion with the values of a, rg, s 
from table 2 and the dotted line the effective range expansion with only a, ro- A) shows the real 
parts and B) the imaginary parts. All amplitudes are in fm. 
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